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Mellin-Barnes (MB) representations have become a widely used tool for the evaluation of Feynman 
loop integrals appearing in perturbative calculations of quantum field theory. Some of the MB 
integrals may be solved analytically in closed form with the help of the two Barnes lemmas which 
have been known in mathematics already for one century. The original proofs of these lemmas solve 
the integrals by taking infinite series of residues and summing these up via hypergeometric functions. 
This paper presents new, elegant proofs for the Barnes lemmas which only rely on the well-known 
basic identity of MB representations, avoiding any series summations. They are particularly useful 
for presenting and proving the Barnes lemmas to students of quantum field theory without requiring 
knowledge on hypergeometric functions. The paper also introduces and proves an additional lemma 
for a MB integral f dz involving a phase factor exp(±i7rz). 
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I. INTRODUCTION 



When working with generalized hypergeometric func- 
tions and relating them through complex contour inte- 
grals, about one century ago, Barnes introduced the two 
following identities: The first Barnes lemma [1], 
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and the second Barnes lemma [2], 
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These identities involve so-called Barnes contour inte- 
grals whose paths of integration run parallel to the imag- 
inary axis and are curved where necessary to ensure that 
they separate positive and negative sequences of poles 
in the integrands. Explicitly, the right poles of gamma 
functions of the type r(a i — z) (at z — ai^ui + 1, . . .) 
must lie to the right of the integration contour, whereas 
the left poles of gamma functions T(ft + z) (at z — 
—ft , —ft — 1, . . .) lie to the left of it. The two Barnes 
lemmas are valid whenever such an integration contour 
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can be found, i.e. unless any of the (Xi + ft is zero or a 
negative integer. Due to the asymptotic behavior of the 
gamma function (see, e.g., Eq. 8.328 1. in [3]), 



\T(x + iy)\ 



\y\- 
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the integrands in (1) and (2) are exponentially damped 
for z — > ±ioo and the integrals converge. 

In the context of loop calculations in quantum field 
theory, Barnes contour integrals arise when using Mellin - 
Barnes (MB) representations by applying the basic iden- 
tity 
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dz B z 



T(-z)T(X + z), (4) 



(A + B)^ T(A) ft loo 2m A*+* 

with the above-mentioned definition of the integration 
contour. In particular, if Re A > 0, the contour can be 
chosen as a straight line in the strip — Re A < Re z < 0. 
When the MB integration in (4) is interchanged with 
other integrations over parameters contained in A or B, 
the convergence of these integrations will usually be en- 
coded in the arguments of the resulting gamma functions. 
The change in the order of integrations is justified if the 
real parts of the arguments of all gamma functions re- 
main positive along the path of the MB integration. Such 
MB integrals may be analytically continued to parameter 
values Re A < 0, provided that the MB integration con- 
tours are deformed in accordance with the requirements 
for Barnes contour integrals described above, separating 
right and left poles from all gamma functions in the in- 
tegrand. 

The identity (4) is valid unless A is a nonpositive in- 
teger. Also the ratio B/A must be away from the nega- 
tive real axis, otherwise the exponential damping of the 
gamma functions for z — > ioo or z — > — ioo is compen- 



sated by (B/A) 2 



and the integral (4) might di- 



verge at one of its boundaries. 
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For \A\ > \B\, the MB integration contour can be 
closed to the right at Rez = +00, and the series of 
residues at z — 0, 1, . . . reproduces the Taylor expan- 
sion of the left-hand side of (4) for |A| > On the 
other hand, for \A\ < \B\, the contour may be closed 
to the left at Rez = —00, and the series of residues 
at z = —A,— A — 1,... coincides with the expansion for 
\A\ < \B\. Therefore the MB representation (4) provides 
a way of separating a sum of two terms A and B without 
requiring a particular hierarchy between them. 

The MB representation (4) was applied originally [4] 
in order to turn massive Feynman propagators l/(k 2 — 
m 2 + iO) into massless ones by separating the mass term 
from the squared momentum, using (4) with A = k 2 and 
B = —m 2 + iO (see also, e.g., [5-7] as early references). 
Later on parametric integrals derived from Feynman di- 
agrams have been evaluated by applying as many MB 
representations as necessary to obtain results in terms of 
gamma functions under the MB integrals (see, e.g., [8- 
11]). The construction of such MB representations for 
Feynman integrals with planar topologies is automated 
by the Mathematica tool AMBRE [12, 13]. 

When evaluating the MB integrations, one has to take 
care of singularities occurring in certain parameter limits 
(like d — ¥ 4 in dimensional regularization). These singu- 
larities arise when a right pole comes arbitrarily close to 
a left pole, pinching the integration contour which sep- 
arates them. Two strategies [14, 15] were formulated to 
systematically extract such singularities analytically. By 
taking a finite number of residues, they both end up with 
integration contours which are not pinched any more be- 
tween right and left poles. These contours can then be 
chosen as straight lines parallel to the imaginary axis 
and the MB integrations may be performed numerically. 
The strategy of [15] is automated in the private computer 
code [16] and by the Mathematica package MB [17]. More 
recently also the strategy of [14] has been implemented 
in the computer program MBresolve [18]. 

MB representations are also employed for the asymp- 
totic expansion of Feynman integrals (see, e.g., [8-11, 19- 
22]). If the MB integral Jdz involves the factor t z , its 
asymptotic expansion in the limit t — > is obtained by 
picking up, order by order, the residues on the right-hand 
side of the integration contour. Similarly the asymp- 
totic expansion in the limit t — > 00 results from the 
residues on the left-hand side of the contour. Czakon's 
code MBasymptotics [23] performs this task automati- 
cally. The contributions obtained from asymptotically 
expanding MB representations often correspond to those 
arising in the strategy of expansion by regions [24-26]. 
This correspondence is emphasized and used for cross- 
checks in [21, 22]. Also combinations of expansion by 
regions with further MB representations for evaluating 
the contributions have proven fruitful [27]. 

The application of MB representations is explained in 
many details in the book [28] . Computer codes are sum- 
marized on the web site [23] . 

When dealing with multifold MB representations, it 



is always preferable to eliminate as many MB integrals 
as possible using the Barnes lemmas (1) and (2) before 
proceeding to the extraction of singularities and numer- 
ical integrations. So these two lemmas have become 
widely used in the application of MB representations. 
They are implemented, in addition to previously men- 
tioned tools, in Kosower's Mathematica package Barnes 
Routines [23]. Let us now turn to the proofs of the two 
Barnes lemmas. 

In his original papers, Barnes focuses on the treatment 
of hypergeometric functions. So, naturally, his proofs 
of the two lemmas employ properties of such functions. 
For proving the first lemma (1), Barnes [1] closes the 
integration contour to one side, which is only possible for 
Re(ai+a2 + /3i + /32) < 1, and sums the two infinite scries 
of residues from the corresponding poles. These series 
are identified with hypergeometric 2-P1 functions of unit 
argument. After relating the 2 -Pi functions to gamma 
functions, using several times T(x) T(l — x) = ir/ sin(7ra;), 
and employing a trigonometric identity, Barnes arrives 
at the right-hand side of (1). The validity of the lemma 
for general complex values of and f$\p, is argued by 
analytic continuation. 

In his proof of the second lemma (2), Barnes [2] starts 
with the series representation of a hypergeometric 3F2 
function of unit argument. The number of gamma func- 
tions depending on the summation index is reduced by 
two using the first Barnes lemma (1) from right to left. 
The series is then identified with a 2-F1 function of unit 
argument which is written in terms of gamma func- 
tions. The resulting expression involves a contour inte- 
gral which can be brought into agreement with the left- 
hand side of (2). This relates the integral from the second 
lemma to the 3F2 function, which, for the specific choice 
of parameters present in (2), however, reduces to a 2-F1 
function of unit argument and can therefore be expressed 
through gamma functions, producing the right-hand side 
of the second lemma. 

These two original proofs by Barnes have been re- 
produced in several books, e.g., [29-31]. Their under- 
standing requires knowledge on the series representa- 
tions and other properties of hypergeometric functions 
which appear in intermediary steps. In the current pa- 
per, however, I present new, elegant proofs for the two 
Barnes lemmas which only use integral transformations 
in a straightforward way, avoiding completely the series 
summations which the original proofs employ. My new 
proofs merely rely on the basic identity of MB repre- 
sentations (4) and on the integral representations of the 
gamma function and Euler's beta function. The latter is 
well known in the field of loop calculations. Therefore, 
when introducing MB representations and the Barnes 
lemmas to students of quantum field theory, I suggest 
to use the proofs presented here rather than the original 
ones by Barnes. 

The new proofs for the two Barnes lemmas are found 
in Sees. II and III, respectively. Sec. IV introduces and 
proves an additional lemma for a MB integral involving 
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a phase factor e ±l ^ z , which, apart from its application 
to loop calculations, is of particular interest for study- 
ing the convergence behavior of such integrals at their 
boundaries ±ioo. Finally conclusions are presented in 
Sec. V. 



II. PROOF OF THE FIRST BARNES LEMMA 

We start with the left-hand side of the first Barnes 
lemma (1), 



dz 
2in 



r(ai - z) T{a 2 - z) r(ft + z) r(ft + z) . (5) 



Let us assume that 

Re(a z + /3j) > Vz, j 



(6) 



such that the integration contour can be chosen as a 
straight line parallel to the imaginary axis with 



Re ft < Rez < Re a, Vi, j . 



(7) 



We insert for T(ai — z) and r(ft + z) the corresponding 
integral representations of the gamma function (see, e.g., 
Eq. 8.310 1. in [3]), 

T(a)= dtt^e' 1 [Rea>0], (8) 
Jo 



and obtain 



dMi e _tl / dt 2 t£ 1 - L e- t * 
(i Jo 

' loo 2^{rj r («2-z)r(/3 2 + z) (9) 



from (5). The change in the order of the integrations is 
justified by the condition (7) for the contour. The basic 
MB identity (4) can easily be reformulated by a shift in 
the integration variable: 
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Inserting (10) for the contour integral in (9) yields 



(10) 



r(« 2 + ft)/ dtj dt 2 ^ ^ 
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o (ii +t 2 ) Q2+fe 
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This double integral factorizes by performing the variable 
transformation t\ = rj£ and t 2 = 7y (1 — £): 



T(a 2 +/3 2 ) / dr?r? 



Ql+^l-l -J) 



X / d^ 
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(i-0 
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While the 77-integral is given by (8), the ^-integral is a 
representation of Euler's beta function (see, e.g., Eqs. 
8.380 1. and 8.384 1. in [3]), 

"W 1 " 1 (1 - = T ^ l)T } a2 ] [Rea t > 0] , 

3 ]"(«!+ a 2 ) 

which is an identity often used in loop calculations. This 
leads to the result 



r(ai + ft) r(ai + f3 2 ) T(a 2 + gi) T(a 2 + /3 2 ) 
r(ai + a 2 + ft + ft) 



(14) 



in agreement with the right-hand side of (1), and thus 
proves the first Barnes lemma. 

The transformations in this proof are valid if the con- 
dition (6) is fulfilled, which is when the integration con- 
tour can be chosen as a straight line (7) or, equivalcntly, 
when the real parts of the arguments of all gamma func- 
tions in the numerator of the result (14) are positive. 
Through analytic continuation, the validity of the first 
Barnes lemma may be extended to general complex val- 
ues of ai i2 and ft j2 , with the exception of parameter 
choices where no integration contour separating right and 
left poles is found in (5), which corresponds to the right- 
hand side (14) being ill-defined as well when at least one 
of the «j -I- ft is a nonpositive integer. 



III. PROOF OF THE SECOND BARNES 
LEMMA 

The left-hand side of the second Barnes lemma (2) 
reads 



3 dz_ r(ai - z) T(a 2 - z) T(ft + z) T(ft + z) 
Q 2iir r(ai + a 2 + ft + ft + ft + z) 
x T(ft + z) . (15) 



As for the other lemma, we assume that the condition (6) 
holds for all pairs of parameters such that the integration 
contour can be chosen as a straight line in the strip de- 
fined by (7). Let us employ the first Barnes lemma (1) to 
replace three of the gamma functions in the integrand: 



r(ft + z)r(ft + z) 

r(ai + a 2 + ft + ft + ft + z) 

1 



r(c*i + a 2 + ft + ft) r(ai + a 2 + ft + ft) 
r 00 ds 

x / — Y{z - s) r(ai + a 2 + ft - s) 

x r(ft + s )r(ft + s ). (16) 

Also this second contour integral over s can be chosen 
along a straight line parallel to the imaginary axis with 

- Reft < Res < Rez < ReaiVi, j . (17) 
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Then the order of the integrations may safely be inter- 
changed and the z-integral yields 



r(c*i - z) r(« 2 - z) r(p! + z) r(-s + z ) 



3 2iir 

r(ai + gi) I> 2 + gi) r(ai - a) r(a 2 - s) 
r(ai + a 2 + /3i — s) 



(18) 



using the first Barnes lemma (1) again. Also the final 
s-integral is solved via the first lemma, 
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2iTT 
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We obtain the result 

r(ai+/3i)r(ai+ / 8 2 )r(ai+/3 3 ) 
T(ai + a 2 + /?i + /3 2 ) r(ai + a 2 + /?i + /3 3 ) 

x r(a 2 + A)r(a 2 + /? 2 )r(Q 2 + /3 3 ) 

r(ai + a 2 + J 2 + 



(20) 



which coincides with the right-hand side of the second 
Barnes lemma (2). So we have proven the second Barnes 
lemma simply by applying three times the first Barnes 
lemma. 

As for the first lemma, the condition (6) can be re- 
laxed to general complex values of the parameters, as 
long as none of the ot^ + f3j is a nonpositive integer. And, 
as before, the arguments of the gamma functions in the 
numerator of (20) reflect the singularity structure of the 
integral (15). 



IV. ADDITIONAL LEMMA FOR 
MELLIN— BARNES INTEGRALS INVOLVING A 
PHASE FACTOR 

As discussed in the introduction (Sec. I), the exponen- 
tial damping of the gamma function towards ±ioo (3) 
usually makes Barnes contour integrals converge well and 
allows for numerical integrations of such integrals. This 
nice feature is potentially spoiled by phase factors e ±l7rz 
in the integrand. Take, for instance, the basic MB iden- 
tity (4) with A = 1 and B = -p + iO (p > 0), thus 
B z = p z e" z . The integrand still vanishes exponentially 
for z — > ioo, but at the lower boundary it behaves as 



P' 



e" 2 r(-z)T(A + z) 



2np 



Re 2 



Im z 



Ro A-l 



(21) 



so the integral only converges for Re A < 0. This di- 
vergence is reflected in the left-hand side of (4), here 
1/(1 — p + i0) A , which exhibits a singular cancelation for 
p ->■ 1 if Re A > 0. 



When an "ordinary" MB integral (which converges at 
z = ±ioo) is solved in closed form, as for the two Barnes 
lemmas (1) and (2), its singularities are encoded in the 
arguments of the resulting gamma functions. These sin- 
gularities originate from right and left poles pinching the 
integration contour between them (see Sec. I). It is de- 
sirable to establish solutions in closed form also for MB 
integrals with phase factors e ±mz where a potential di- 
vergence at z = ±ioo is parametrized by the resulting 
gamma functions as well. 

One such identity is provided by the following addi- 
tional lemma: 



dz e±MZ r(a - z) Tijii + z) rCgg + z) 
2i7T r(7 + z) 



, i7ra r(a + /? 1 )r(a + ( s 2 )r( 7 . 
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r(7-ft)r( 7 -0 2 ) 



(22) 



While the first two gamma functions in the numerator of 
the right-hand side of (22) reflect the singularities origi- 
nating from right and left poles, the third gamma func- 
tion parametrizes the divergence at z = =pioo. In the case 
of the phase factor e +t7TZ , the integrand in (22) exhibits 
the asymptotic behavior 



T{a-z)Y{(3 l +z)T{p 2 + z) 



T( 7 + z) 

2 7 r|Imz| Rc ( Q+/3l+ft -^- 1 . (23) 
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So the integral only converges at z — —zoo for 

Refr-a-ft-ft) >0, (24) 

which matches the argument of the corresponding gam- 
ma function in the right-hand side of (22). 

Let us turn to the proof of the additional lemma (22). 
In analogy to the previous proofs we assume 



Re(a + /3j) > OVj 



(25) 



such that the integration contour can be chosen as a 
straight line parallel to the imaginary axis with 



Re/3j < Rez < ReaVj . 



(26) 



Here we also have to require the condition (24) for the 
convergence of the integral at both boundaries z — ±ioo. 
We start with the left-hand side of (22) and replace two 
gamma functions using (13): 



1 



r(fe + z) = 
r( 7 + Z ) r( 7 - (3 2 ) J 



02+Z-l 



(i-O 



(27) 



The convergence of the integral (27) follows from (24) 
and (26). We may safely change the order of integrations 
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and perform the contour integral first employing the basic 
MB identity in the form (10): 



d^ 
2i?r 



C T(a- z)T{f3 1 +z) 



(-£±i0)« 

T(a + /3i) (-e±io) a (i-o _a_A 



Then the ^-integral is given by (13): 



dee 



(i-0 



7-a-0i-/3 2 -l 



r(7-j9i) 



(28) 



(29) 



By assembling the pieces from (27), (28), and (29) we 
obtain the right-hand side of (22). 

The validity of the additional lemma (22) can be ex- 
tended via analytic continuation to general complex val- 
ues of a, /3i,2> and 7 by relaxing the condition (25), as 
discussed for the two previous lemmas. However, the con- 
dition (24) is crucial for the convergence of the contour 
integral and cannot be waived. 



V. CONCLUSIONS 

New, elegant proofs have been presented for the two 
Barnes lemmas which are often used in loop calculations 
in order to solve MB integrals. The new proofs avoid the 
series summations and representations of hypergeometric 
functions which are used in the original proofs by Barnes. 
They are well suited for courses in quantum field theory 
treating MB representations because they are only based 
on identities generally known there. 

Also an additional lemma has been introduced and 
proven which parametrizes the divergence of a MB in- 
tegral J dz when a phase factor e ±l7TZ is present. This 
new lemma has been used by myself in many evaluations 
of Feynman integrals, and other authors have probably 
solved this and similar MB integrals as well, e.g. by taking 
a series of residues and summing them up via a hyperge- 
ometric function. But, to my knowledge, this lemma has 
not been published in closed form so far. 
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